clear;
m = [5:2:50];
r = [1,5:5:200];
lambda = 700e-9;

[xa,yr] = meshgrid(m,r);
Rd = zeros(size(xa));
%%
for i = 1:size(Rd,1)
    for j = 1:size(Rd,2)
%        Rd(i,j) = Reflectance_Diffuse(yr(i,j),xa(i,j),lambda);
        Rd(i,j) = Reflectance(yr(i,j)*1e-9,xa(i,j),lambda);
    end
end
%%
mesh(xa,yr,Rd);
set(gca,'FontSize',20);
xlabel('m_{\Delta}','FontSize',20);
ylabel('r_{\Delta}(nm)','FontSize',20);
zlabel('R/R_0','FontSize',20);